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RESUMEN 

Anunciamos la publicacion de un catalogo de modelos fisicos de discos 
de acrecion irradiados por la estrella central alrededor de estrellas jovcncs, 
basados en las tecnicas de modelizacion de D'Alessio ct al. El catalogo WWW 
incluye ~ 3000 modelos de disco para diferentes estrellas centrales, tamafios 
del disco, inclinaciones del disco, contcnidos de polvo y tasas de acrecion de 
masa. Para cada uno de ellos, los perfiles radiales de las propicdades fi'sicas 
del disco y distribuciones espectrales de energfa sinteticas se pueden consul- 
tar o descargar desde la pagina web para compararlas con observaciones. El 
catalogo est a accesible en [ht tp : / / www-cf a . harvard . edu/ youn gstars/ dalessio/| 
(EEUU) , |ht tp : / / www . astrosmo . unam . mx/ ^dalessio/ (Mexico), y en 
[http : //www ■ laef f . esa . es/mo dels/dalessio/ (Espana). 

ABSTRACT 

We announce the release of a catalog of physical models of irradiated 
accretion disks around young stars based on the modelling techniques by D'Alessio 
et al. The WWW catalog includes ~ 3000 disk models for different central 
stars, disk sizes, inclinations, dust contents and mass accretion rates. For any of 
them, radial profiles of disk physical parameters and synthetic spectral energy 
distributions can be browsed and downloaded to compare with observations. 
It can be accessed at http : //www-cf a .harvard. edu/youngst ars/dalessio7| 
(US), l^ tp : //www, astrosmo . unam.mx/ ^dalessio/| (Mexico), and at 
[http: //www . laef f . esa.es/models/dalessio/ (Spain). 
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1. INTRODUCTION 

The old idea that stars arc born surrounded by 
disks, which may form planetary systems, has found 
strong observational support in the last couple of 
decades. The properties of these disks arc quantified 
from the comparison between different observations 
and models. For instance, disk mass accretion rates 
have been inferred from the analysis of the short 
wavelength excess, modeled as produced by accre- 
tion shocks at the stellar surface (Hartigan, Edwards 
& Ghandour 1995; Gullbring et al. 1998; Calvet & 
Gullbring 1998; Hartmann ct al. 1998). Also, mass 
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accretion rates and inner disk radii have been esti- 
mated from models of different line emission profiles 
thought to form in the magnetospheric flows con- 
necting the disks to their central stars (MuzcroUc ct 
al. 1998a,b and c; MuzcroUc ct al. 2001). Kinetic 
information on disks, central star masses, details of 
the vertical temperature distribution and molecular 
abundances are quantified from detailed analysis of 
different molecular lines (e.g., Dutrey et al. 1996; 
Dutrey, Guilloteau & Guehn 1997; Guilloteau & 
Dutrey 1998; Simon, Dutrey & Guilloteau 2000; Na- 
jita et al. 2000; Dartois, Dutrey & Guilloteau 2003; 
Aikawa et al. 2002, 2003; Qi et al. 2003; Carr, Toku- 
naga & Najita 2004). Other disk properties, such as 
masses, degree of flaring of the disk surface and dust 
properties, are inferred from the spectral energy dis- 
tributions (SEDs) of young stars, from near IR to 
radio-frequencies, using different kinds of disk mod- 
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els (e.g., Kenyon & Hartmanii ll987t Beckwith et al. 
1990; Calvet et al. 1991, 1992; Malbet & Bertout 
1991; Beckwith & Sargent 1991; Chiang & Goldre- 
ich 1997, 1999; D'Alessio et al. 1998, 1999, 2001; 
DuUemond, Dominik & Natta (2001); DuUcmond, 
van Zadelhoff & Natta (2002); DuUemond (2002); 
Malbet, Lachaume & Monin 2001; Lachaume, Mal- 
bet & Monin 2003). These models go from sim- 
ple power-law descriptions of the disk mass surface 
density, temperature and opacity, to detailed mod- 
els where different heating mechanisms are included, 
and where the radiative transfer and dust opacity 
are calculated with different degrees of sophistica- 
tion. Disks have been also imaged at near IR (e.g., 
Stapclfcldt et al. 1998; Koresko 1998; Malbet et al. 
1998; Weinberger ct al. 1999; Padgett et al. 1999; 
Tuthill et al. 2002; Colavita et al. 2003), mid IR 
(e.g., McCabe, Duchene, & Ghez 2003) and millime- 
ter wavelengths (Dutrey et al. 1996; Wilner, Ho & 
Rodriguez 1996; Rodriguez et al. 1998; Mannings & 
Sargent 2000). Analysis of these observations with 
different models has given information about the ra- 
dial and vertical structures of the disks. In general, 
models with some degree of complexity are required 
to fit multi-wavelength observations of a given ob- 
ject (e.g.. Lay, Carlstrom & Hills 1997; Wilner & 
Lay 2000). 

In the present paper we describe a database with 
a series of models of disks around pre-main sequence 
stars. These models, which are described in detail in 
Merfn (2004), are self-consistently calculated given 
the properties of the central star, the disk and its 
dust, usin g the m eth ods d eveloped by D'Alessio et 
al. l|l998L Il999l and l200lj) . We make these mod- 
els available with the hope that they can be use- 
ful for fitting observations of young stars, helping to 
extract more information than simpler and usually 
faster approaches. In a forthcoming paper (Merfn 
et al. 2005), relationships between disk, star and 
dust parameters and observable quantities such as 
spectral indices, colors, etc., will be presented and 
discussed in detail, to enable a faster search for the 
best model fit to the observations of a given disk or 
a survey of disks. In addition, comparison of models 
in this database with other modeling efforts will be 
helpful to test the underlying assumptions in each set 
of models. Finally, we hope that this database will 
help increase our knowledge of the intrinsic proper- 
ties of the disks around young stars. 

2. DISK MODELS 

The di sk models from D'Alessio et al. (|l998lll999l 
and l200l|) have the following characteristics: 



• Energy is transported by radiation, convection 
and a turbulent flux, the first one being the most 
important mechanism given the disk parameters 
we have explored in this work. 

• The disk mass surface density distribution is a 
consequence of conservation of angular momen- 
tum flux given the functional form of the viscos- 
ity coefficient and the disk mass accretion rate. 

• The gas and dust are in thermal balance, and 
are heated by viscous dissipation, stellar irradia- 
tion, and ionization by energetic particles (from 
cosmic rays and radioactive decay). For the set 
of parameters of models in the database, vis- 
cous dissipation is important in regions close to 
the star and at the disk midplane, while stellar 
irradiation is important for most of the disk. 

• The fraction of the stellar radiative flux in- 
tercepted by the disk is evaluated at the sur- 
face where the mean radial optical depth to the 
stellar radiation is unity. This surface is self- 
consistently calculated given the disk structure 
and properties of its dust, and the intercepted 
flux is used as a boundary condition in the in- 
tegration of the radiative transfer equation. 

The models are based on the following simplify- 
ing assumptions: 

• The disk is in steady state, with a mass accre- 
tion rate A/= dM/dt. 

• The disk is geometrically thin: H/R << 1, 
where H is the gas scale height of the disk and 
R is the radial distance. With this assumption, 
the vertical and radial structures can be calcu- 
lated separately. 

• The viscosity coefficient is given by v = a H Cg, 
following the g -prescription from Shakura & 
Sunyaev lll973l) . where is the local sound 
speed and a is the viscosity parameter, assumed 
to be constant through the disk. The turbu- 
lent flux of energy is calculated consistently with 
the a prescription, assuming a Prandtl number 

Pr = l. 

• Dust and gas are well mixed in the entire disk, 
i.e., the dust to gas mass ratio of each dust ingre- 
dient and the grains size distribution are both 
taken to be constant in those regions where the 
temperature is lower than the sublimation tem- 
perature. 
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• Dust and gas are in thermal equilibrium and a 
unique temperature is calculated for both com- 
ponents. 

• The radiation field is considered in two sepa- 
rate regimes: stellar radiation (UV, optical and 
near IR) and disk radiation (IR to radio wave- 
lengths), and we use mean opacities calculated 
with appropriate weighting functions to describe 
the interaction of both fields with the disk ma- 
terial. 

• The radiative transfer in each regime is done by 
solving the two first moments of the radiative 
transfer equation. 

• The inner disk is truncated at the dust de- 
struction radius and the inner cylindrical sur- 
face ("the wall") emits as a blackbody of 140 
K (Natta et al. 2001, MuzeroUe et al. \200^ . 
This is a new feature added to the models, which 
has not been included in the models described in 
the previous papers (D'Alessio et al. 1998, 1999, 
2001). We follow the prescriptions in D'Alessio 
et al. (2003) to calculate the emission from the 
wall, taking into account the contributions of 
the stellar and accretion luminosities to calcu- 
late the position of the dust destruction radius, 
and the dependence on inclination of the emit- 
ting area. More specifically, we compute the in- 
ner disk radii with eq. (3) from MuzeroUe et al. 
(2003) where we take qxd/ i^d = 2.5, which gives 
the best fits to the excess near-IR continuum 
in Classical T Tauri disks. Then, we assume a 
fixed ratio between the height of the wall and 
the scale height at the dust destruction radius, 
Zwaii = 4i?. However, we list separately the disk 
and wall contributions, allowing the user to in- 
clude different heights or radii for the wall. 

2.1. Comparison with observations 

The synthetic SED of a model calculated for typ- 
ical stellar parameters, mass accretion rate, disk ra- 
dius, inclination angle and a dust grain size distri- 
bution with millimeter-sized grains, fits the median 
observed SED of the Classical T Tauri Stars in th e 
Taurus Molecular Cloud (see D'Alessio et al. 1200 1|) . 
Also, Merm et al. (2004a) used the models to fit 
the detailed SEDs of two HAeBe stars (namely HD 
34282 and HD 141569) and showed that not only the 
age was responsible for the disk evolution but also 
that the metallicity could play an important role. 

Allen et al. (2004) have constructed color-color 
diagrams for different clusters using observations 



from IRAC-Spitzer. They find three distinctive re- 
gions in the diagrams: one corresponding to stars 
without IR excess, other for Classical T Tauri Stars 
and another for embedded objects. The synthetic 
IRAC colors of a subset of models from this library 
(T* = 4000 K and age=l Myr) agree very well with 
those of observed Classical T Tauri Stars (see also 
Siciha- Aguilar et al. 2004; Hartmann et al. 2004). 

It is important to mention that it is very difficult, 
if not impossible, to constrain the model parame- 
ters by considering only the SED of a given object. 
The main problem is that the SED is given by the 
monochromatic flux emerging from whole disk, i.e., 
it is a spatially integrated (luantity. Thus, detailed in- 
formation on the spatial distribution of the emergent 
intensity, which is more directly related to the disk 
structure, is hidden in the fiux. Fortunately, different 
parts of the SED are sensitive to different combina- 
tion of disk parameters, but more than one of these 
combinations can produce the same emergent flux 
from the disk in a particular wavelength range. For 
instance, the near IR SED is dominated by the emis- 
sion of the wall at the dust sublimation radius. For 
a given sublimation temperature, different silicate 
composition, grain sizes, total (stellar and accretion 
shock) luminosities result in different values of the 
sublimation radius, but different inclination angles 
and wall heights can be assumed to obtain a similar 
wall SED. In principle, a SED which spans a wide 
range of wavelengths (from UV to radio-frequencies) 
might help to disentangle some of the disk proper- 
ties when a self-consistent model is used. For in- 
stance, the continuum at mid-IR emerges from opti- 
cally thick zones of the disk. In this spectral range, 
the emergent flux depends on the disk mass accre- 
tion rate (if the accretion luminosity is similar or 
larger than the stellar luminosity), the disk incli- 
nation angle, and the grains composition and size 
distribution (which determines the height of the ir- 
radiation surface). The 10 and 18 /zm silicate bands 
reflect the silicate composition and grain sizes in the 
atmosphere, thus a self-consistent model should have 
the same type of grains producing the observed fea- 
tures and absorbing the stellar radiation which de- 
termines the disk vertical temperature distribution 
and its mid-IR emergent continuum. However, not 
all the degeneracy can be removed by modeling the 
SED alone, specially if the dust at the midplane in 
the outer disk has a different size distribution than 
the dust in the atmosphere in the inner disk, which 
is expected if dust is settling and growing in disks. 
The only way to remove all possible degeneracies is 
when the SED and images if a self-consistent model 
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are compared to multi-frequency high angular reso- 
lution observations. In this first release of the catalog 
we are making available only SEDs, but intensity dis- 
tributions at different wavelengths can be calculated 
upon request for a given set of model parameters. 

3. WEB-BASED MODEL LIBRARY 

We have constructed a grid of models for different 
central stars and a range of values for the physical 
parameters of the disk and its dust. 

For the central stars, we chose 17 stars with spec- 
tral types from K7 to B9 {T^s = 4000-10000 K re- 
spectively) and ages of 1 and 10 Myrs. Stellar param- 
eters were take n from the pre-main sequence tracks 
by Siess et al. (|2000() . 

For the disk, we considered four val- 
ues of the accretion rate, namely M = 
10~^ 10"^ 10"^ lO-^Mg/yr, a typical viscos- 
ity parameter a = 0.01, three values of the disk 
radius i?disk = 100, 300, 800 AU and two inclination 
angles i = 30, 60 degrees. 

For the dust we use t he abundances proposed 
by Pollack et al. a grain size distribu- 

tion given by n{a) = a~P with p = 3.5, 2.5, with a 
minimum size typical of interstellar grains (a„ii,i = 
0.005 ^m) and six different maximum grain sizes 
amax = l, 10, 100 /im, 1mm, 1cm, 10 cm. 

Columns (1) to (6) of table ^ show the stellar 
effective temperatures, spectral types, ages, radii, 
masses and luminosities for the central stars in the 
grid of disk models, column (7) lists the mass ac- 
cretion rates of the disks, column (8) has the ac- 
cretion luminosity, and columns (9) and (10) show 
the wall radius, in astronomical units and in stellar 
radii. Notice that, for a given central star, the wall 
radius increases with mass accretion rate, reflecting 
the fact that we are including the irradiation from 
the accretion shocks at the stellar surface as an ad- 
ditional heating source of the dust (see MuzeroUe et 
al. 2003). 

The disk models were computed in parallel with 
Linux PCs at the Harvard-Smithsonian CfA (MA, 
US), LAEFF (Madrid, Spain) and the Linux cluster 
nostromo at the CRyA (Morelia, Mexico) and used 
a total CPU time of approximately 8000 hours. 

4. CONCLUDING REMARKS 

The WWW catalog contains currently more than 
3000 model SEDs and disk structures. We hope this 
library of models will be a valuable tool for analyz- 
ing the SEDs of young stars surrounded by accretion 
disks. It is a dynamical website. We plan to continue 
filling the catalog with the emission maps and new 



disk models with different degrees of dust settling to 
the midplanc. 
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TABLE 1 

PARAMETERS FOR THE CENTRAL STARS AND DISKS IN THE LIBRARY OF DISK MODELS. 



Central Star 



Irradiated accretion disk model 





SoTvDe Aee 


L, 






M 




I^wall 




K 


Myr 




-LV,0 


XVJ.0 


Mr-, vr""' 




AU 


R* 


4000 


K7 1 


1.60 


2.64 


0.70 


1 X 10"^ 


0.009 


0.11 


8.65 












1 X 10~* 


0.087 


0.11 


8.88 












1 X 10"^ 


0.869 


0.13 


10.68 












1 X 10^^ 


8.694 


0.27 


21.58 


4000 


K7 10 


0.31 


1.16 


0.80 


1 X 10"'' 


0.022 


0.05 


8.97 












1 X 10~* 


0.217 


0.06 


11.31 












1 X 10"^ 


2.167 


0.13 


24.52 












1 X 10"'' 


21.67 


0.39 


73.05 


4500 


K4 1 


3.76 


3.20 


1.40 


1 X 10"'"* 


0.014 


0.16 


10.92 












1 X 10~* 


0.137 


0.17 


11.15 












1 X 10~^ 


1.375 


0.19 


12.80 












1 X 10^*^ 


13.75 


0.35 


23.63 


4500 


K4 3 


1.60 


2.08 


1.34 


1 X 10~^ 


0.020 


0.11 


11.06 












1 X 10~* 


0.202 


0.11 


11.67 












1 X 10"^ 


2.024 


0.16 


16.54 












1 X 10"*^ 


20.24 


0.39 


40.61 


4500 


K4 10 


0.60 


1.28 


1.10 


1 X 10~^ 


0.027 


0.07 


11.18 












1 X 10~* 


0.270 


0.08 


13.17 












1 X 10~^ 


2.700 


0.15 


25.65 












1 X 10"'" 


27.00 


0.44 


74.19 


5000 


Kl 1 


14.86 


5.14 


3.00 


1 X 10"^ 


0.018 


0.32 


13.56 












1 X 10^* 


0.184 


0.33 


13.64 












1 X 10^^ 


1.834 


0.34 


14.37 












1 X 10"*^ 


18.34 


0.48 


20.26 


5000 


Kl 10 


1.72 


1.75 


1.40 


1 X lO"'' 


0.025 


0.11 


13.64 












1 X 10"** 


0.251 


0.12 


14.50 












1 X 10"^ 


2.514 


0.17 


21.25 












1 X 10"^ 


25.14 


0.43 


53.52 


6000 


GO 1 


59.10 


7.13 


3.50 


1 X lO"'' 


0.015 


0.65 


19.49 












1 X 10"* 


0.154 


0.65 


19.51 












1 X 10"^ 


1.542 


0.65 


19.74 












1 X 10"^ 


15.42 


0.73 


21.88 


6000 


GO 10 


5.91 


2.25 


1.60 


1 X lO"'' 


0.020 


0.20 


19.56 












1 X 10"* 


0.223 


0.21 


19.89 












1 X 10"^ 


2.234 


0.24 


22.92 












1 X 10"^ 


22.34 


0.45 


42.70 



6 



D'ALESSIO ET AL. 



TABLE 1 
(CONTINUED). 



Central Star 



T^ff 

cri 


SoTvDe 




L, 




M* 


M 




I^wall 


r^wall 


K 




Mvr 










-LJ0 


AU 


R, 


7000 


Fl 




130.39 


7.78 


4.00 


1 X 10~^ 


0.016 


0.96 


26.53 














1 X 10~* 


0.161 


0.96 


26.54 














1 X 10~^ 


1.615 


0.97 


26.69 














1 X 10~® 


16.15 


1.02 


28.12 


7000 


Fl 


10 


11.02 


2.20 


1.70 


1 X 10"^ 


0.024 


0.28 


27.30 














1 X 10^® 


0.243 


0.28 


27.57 














1 X 10~^ 


2.428 


0.31 


30.13 














1 X lO""" 


24.28 


0.50 


48.81 


8000 


A6 


1 


165.08 


6.70 


4.0 


1 X 10^^ 


0.019 


1.08 


34.71 














1 X 10"® 


0.188 


1.08 


34.71 














1 X 10"'' 


1.879 


1.09 


34.90 














1 X lO"*" 


18.79 


1.14 


36.63 


8000 


A6 


10 


12.62 


1.85 


1.9 


1 X 10~^ 


0.032 


0.30 


34.75 














1 X 10~® 


0.322 


0.30 


35.15 














1 X lO"'^ 


3.227 


0.33 


38.89 














1 X 10"'' 


32.27 


0.56 


65.46 


9000 


A2 




215.39 


6.05 


4.0 


1 X lO"'' 


0.021 


1.23 


43.84 














1 X 10~® 


0.207 


1.23 


43.86 














1 X 10"'' 


2.077 


1.24 


44.05 














1 X lO"'' 


20.77 


1.29 


45.91 


9000 


A2 


3 


71.00 


3.47 


2.7 


1 X 10"' 


0.024 


0.71 


43.90 














1 X 10"** 


0.244 


0.71 


43.94 














1 V 1 n^'^ 




n 79 


AA RA 














1 X 10"^ 


24.45 


0.82 


50.89 


9000 


A2 


10 


17.08 


1.70 


2.0 


1 X 10"' 


0.037 


0.34 


43.99 














1 X 10"* 


0.370 


0.35 


44.41 














1 X 10"^ 


3.697 


0.38 


48.46 














1 X 10"*^ 


26.97 


0.62 


78.16 


10000 


B9.5 


1 


251.94 


5.30 


4.0 


1 X 10"' 


0.024 


1.33 


54.13 














1 X 10"** 


0.237 


1.33 


54.15 














1 X 10"^ 


2.371 


1.34 


54.38 














1 X 10"^ 


23.71 


1.40 


56.62 


10000 


B9.5 


10 


29.20 


1.80 


2.3 


1 X 10"' 


0.040 


0.45 


54.30 














1 X 10"® 


0.402 


0.46 


54.63 














1 X 10"^ 


4.015 


0.48 


57.87 














1 X 10"*^ 


40.15 


0.70 


83.62 
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